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The Laser Interferometer Gravitational Wave Observatory (LIGO) has performed a third science 
run with much improved sensitivities of all three interferometers. We present an analysis of ap- 
proximately 200 hours of data acquired during this run, used to search for a stochastic background 
of gravitational radiation. We place upper bounds on the energy density stored as gravitational 
radiation for three different spectral power laws. For the flat spectrum, our limit of fio < 8.4 x 10~* 
in the 69-156 Hz band is ~10^ times lower than the previous result in this frequency range. 

PACS numbers: 04.80.Nn, 04.30.Db, 95.55.Ym, 07.05.Kf, 02.50.Ey, 02.50.Fz, 98.70.Vc 



A stochastic background of gravitational waves could 
result from the random superposition of an ex- 
tremely large number of unresolved and independent 
gravitational-wave (GW) emission events. Such a back- 
ground is analogous to the cosmic microwave background 
radiation (CMBR), though its spectrum is unlikely to be 
thermal. The emission events could be the result of cos- 
mological processes, as with the CMBR, but occurring 
much earlier after the big bang — e.g., during inflation. 
The events could also be due to more recent astrophys- 
ical processes. Placing upper limits on or detecting the 
energy density of a stochastic background of gravitational 
waves is one of the long term goals of GW detectors. 

The stochastic background spectrum is typically char- 
acterized in terms of a dimensionless quantity Q.^{f): 
the GW energy density per unit logarithmic frequency, 
divided by the critical energy density pc to close the 
universe. The critical density, and thus figw(/), de- 
pend on the Hubble expansion rate Hq\ in this letter, all 



bounds on rig„(/) will be for = 72 km sec^^ Mpc^^ 
[ij. We search for power laws of the form ilgw(/) = 
ila{f /lOO Hz)". The choices of a are motivated by po- 
tential source models 0: (a = 0) predicted by inflation- 
ary or cosmic string models; (a = 2) rotating neutron 
stars; {a = 3) pre-big-bang cosmology. 

Previous direct measurements of a stochastic back- 
ground, in the ^10 Hz to ^10'* Hz band accessible to 
earth-based detectors, have been limited to establishing 
upper limits on r2gw(/) greater than unity, with the best 
and most recent result using LIGO's first science data 
finding f^o < 44 At much lower frequencies, space- 
craft Doppler tracking has established JIq < 0.027 in the 
band 10^^ — 10""^ Hz and radio pulsar timing has 
bounded fJo < 10^ in a decade band around 10~* Hz 0. 
In this paper we report new limits on a stochastic GW 
background for frequencies around 100 Hz, using data 
from the LIGO GW detectors. In terms of GW energy 
density, these limits are nearly five orders of magnitude 
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below previous limits in this frequency band. 

LIGO is composed of three GW interferometers at two 
sites: the 4 km HI and 2 km H2 detectors, collocated 
at Hanford, WA; and the 4 km LI detector, located in 
Livingston Parish, LA 0. Each detector is a power- 
recycled Michelson interferometer, with 4 km (or 2 km) 
long Fabry-Perot cavities in each of its orthogonal arms. 
These interferometers are sensitive to quadrupolar oscil- 
lations in the space-time metric due to a passing GW, 
measuring directly the GW strain amplitude. While 
the detectors are still being commissioned to perform at 
their designed sensitivity, several dedicated data collec- 
tion runs have been performed. The detector configu- 
ration and performance during LIGO's first science run 
(SI) is described in Ref. 0- In this analysis we use data 
from the third science run (S3), carried out from 31 Oct 
2003-9 Jan 2004, with significantly improved sensitivity 
compared to previous runs. Fig. ^ shows reference am- 
plitude spectra of equivalent strain noise for S3. These 
were typical noise levels for the detectors, though there 
was also significant variation during the run; e.g., HI 
displayed a general trend of improved performance over 
time, with a minimum noise level at the end of the run 
approximately a factor of 2 below that in Fig. 
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FIG. 1: Reference sensitivity curves during the S3 data run, 
in terms of equivalent strain noise density. Also shown is 
the f~^^^ strain noise level corresponding to the upper limit 
found in this analysis, Qq = 8.4 x 10"*, and the strain noise 
goal for the two 4-km interferometers. 

To search for a stochastic background, we cross- 
correlate the strain data from a pair of detectors, taking 
advantage of the fact that the instrumental noise of one 
interferometer will, in general, be uncorrelated with that 
of the other interferometers. This is more clearly the 
case for the widely separated interferometer pairs (i.e., 
LI paired with HI or H2), for which there are only a 
few paths through which instrumental correlations could 
arise. The H1-H2 cross-correlation, on the other hand, is 
susceptible to a much broader range of potentially corre- 
lated instrumental noise sources. This letter presents an 



analysis of the inter-site pairs only; the H1-H2 measure- 
ment, though off'ering potentially higher sensitivity due 
to their collocation, requires additional techniques to ad- 
dress instrumental noise, and may be presented in a later 
publication. 

The cross-correlation is performed in the frequency do- 
main, using a linear filter that optimizes the expected 
signal-to-noise ratio, given the detectors' noise spectra 
and the targeted stochastic background (see Ref. Q and 
references cited within). Specifically, with si(/) and 
S2(/) representing the Fourier transforms of the strain 
outputs of two detectors, the cross-correlation is com- 
puted as: 



Y = 



df 



df STif~r)§uf)Qif')hin , (1) 



where St is a finite-time approximation to the Dirac delta 
function. The optimal filter Q has the form: 



0(/)=AA 



7(/)^gw(/) 



(2) 



where TV is a normalization factor. Pi and P2 are the 
strain noise power spectra of the two detectors, S'gw is the 
strain power spectrum of the stochastic background being 
searched for (S'gw(/) = {iH"^ /lG-K^)f-^n^^{f)), and the 
factor 7 is called the overlap reduction function . This 
factor, defined so that its absolute value is at most unity, 
gives the frequency variation of the cross-correlation aris- 
ing from an isotropic stochastic background, for sepa- 
rated or non-aligned detectors (7(7) = 1 at all frequen- 
cies for the co-located detector pair, HI and H2). 

The optimal filter is derived assuming that the intrin- 
sic detector noise is Gaussian and stationary over the 
measurement time, uncorrelated between detectors, and 
uncorrelated with and much greater in power than the 
stochastic GW signal. Under these assumptions the ex- 
pected variance, Oy , of the cross-correlation is dominated 
by the noise in the individual detectors, whereas the ex- 
pected value of the cross-correlation Y depends on the 
stochastic background power spectrum: 

rp /"OO 

4 ^ _ ^ / |Q(/)|2P2(/) (3) 



{Y)=T df ^{f)S^^{f)Q{f) , (4) 

JO 

where T is the duration of the measurement. 

Analysis Details.-The analysis is implemented simi- 
larly to the method detailed in Ref. The data set 
from a given interferometer pair is divided into equal- 
length intervals, and the cross-correlation Y and theo- 
retical cry are calculated for each interval, yielding a set 
{Yj,(Jy,} of such values, with / labelling the intervals. 
This data segmentation is useful for dealing with long- 
term non-stationarity of the detector noise, by choosing 
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an interval length over which the noise is relatively sta- 
tionary. The interval length for this analysis is 60 sec. 
The cross-correlation values are combined to produce a fi- 
nal cross-correlation estimator, l^pt, that maximizes the 
signal-to-noise ratio, and has variance (T^p^^'. 

The normalization factor M is defined such that the point 
estimate of J7a and its standard deviation are given by: 

Before computing the cross-correlation, each 60 sec 
data interval is decimated (from 16384 Hz to 1024 Hz), 
high-pass filtered (40 Hz cut-off) and Hann- windowed. 
The windowing step protects against spectral leakage of 
strong lines that may be present in the data, but at the 
same time a Hann window reduces the effective length 
of the interval by nearly a factor of 2 (approximately 
the mean value of the Hann window). To recover the 
loss in signal-to-noise, the data intervals are overlapped 
by 50%, so that each data point receives full weighting 
in the analysis (except for end effects). This introduces 
some correlation between the Yj and ctyj for adjacent 
values of /, complicating the formulae for the optimal 
estimator and its variance, derivations of which may be 
found in Ref. 0- 

The detectors' strain noise power spectral densities 
(PSDs) are estimated for each interval in order to cal- 
culate the optimal filter Q[f). Welch's modified peri- 
odogram method of PSD estimation is used, averaging 
58 periodograms formed from 4 sec-long, 50% overlap- 
ping data windows. The PSD for interval / is formed 
from the two 60 sec data intervals preceding and fol- 
lowing, but excluding the data within, interval /. This 
technique eliminates a bias (underestimate) in the cross- 
correlation that would otherwise exist, due to non-zero 
covariance between the s\-S2 cross-spectrum and the cor- 
responding power spectra. However, short-term changes 
(typically increases) in detector noise may produce out- 
liers, because excess instrumental noise within interval / 
is not reflected in its PSD. This is addressed by apply- 
ing a consistency test on cry^ over 3 consecutive intervals; 
those intervals for which the cry, differ by more than 20% 
are eliminated from the analysis. For this analysis, ap- 
proximately 20% of the data were rejected by this cut, 
leaving 218 hours. 

To compute the cross-correlation (Eq. QJ, the raw 
detector data are calibrated, in the frequency domain, 
into strain units using interferometer response functions. 
These functions are calculated once per 60 sec, using a 
measurement of the response of an interferometer to a 
sinusoidal calibration force, averaged over 60 sec. The 
frequency domain values of {si, S2}, given at a frequency 
spacing of 1/60 Hz, are binned to the resolution of the 
optimal filter (frequency spacing of 1/4 Hz), and the in- 
tegrations in Eqs.^andl^lare performed for the different 
Sgw if) ■ 
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FIG. 2: Coherence between HI and LI during S3, showing 
a few small, but significant, coherent peaks at multiples of 
16 Hz. The horizontal line corresponds to the statistical ex- 
pectation value of l/Na^vg, = 3.3 x 10"^, where A'avg is the 
number of periodogram averages. The inverted triangles at 
the top of the graph indicate the discrete frequencies omitted 
from the analysis. The inset histogram shows that the coher- 
ence values (r) follow the expected exponential distribution 
(dashed line). 



As was important in the earlier analysis of Ref. j^], 
frequency bins corresponding to known or potential in- 
strumental correlation artifacts are excluded from the 
frequency domain integrations. An obvious example of 
inter-site correlations comes from the 60-Hz AC supply 
lines used to power the detectors. The 60-Hz modula- 
tion and its harmonics are present to some degree in the 
detector electronics, and thereby infiltrate the strain out- 
put signal (as Fig.^shows). Between LI and HI, 2, the 
power lines tend to be well correlated over time scales 
shorter than ~100 sec, with ever-decreasing correlation 
over longer times. To exclude the possibility of any resid- 
ual long-term power line correlation, the (60 Hz, 120 Hz, 
. . . ) bins are excluded from the integration. 

Another narrow-band source of instrumental correla- 
tion stems from imperfections in the detectors' data ac- 
quisition systems. The peaks in Fig. 12 at multiples of 
16 Hz, were produced by slight but periodic corruption 
of the data at each site. The data acquisition timing at 
each site is controlled by clocks synchronized to the 1 
pulse-pcr-second signals produced by Global Positioning 
System (GPS) receivers. The 16 Hz periodicity of the 
data corruption was controlled by these clocks, resulting 
in persistent inter-site correlations at multiples of 16 Hz. 
Thus, the (16 Hz, 32 Hz, ... ) bins are also excluded from 
the integration. (After S3, the offending clock modules 
were identified and the problem was corrected.) 

Results.~FoT each power law searched for, we calculate 
the optimal cross-correlation statistic and its variance 
(Eqs.O. The Hl-Ll results are summarized in Table 
The H2-L1 correlation is approximately a factor of 7 less 
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Power 


Frequency 


Pt Estimate 


Statistical Error 


Calibration Error 


Upp 


er Limits 


Law 


Range 






(HI) 


(LI) 


ng„(/) 


5g4 (,/) (Hz ^/") 


a = 


69 - 156 Hz 


-6.0 X 10"'' 


7.0 X IQ-" 


±9% 


±15% 


8.4 X lO""* 


1.2 X 10"2=* (//lOO Hz)"^/2 


a = 2 


73 - 244 Hz 


-4.7 X 10"'* 


7.2 X 10"* 


±9% 


±15% 


9.4 X 10"**(//100 Hz)^ 


1.2 X lO"^'' if /WO Hz)-^/2 


a = 3 


76 - 329 Hz 


-4.0 X 10"" 


6.2 X 10""* 


±9% 


±15% 


8.1 X 10""* if /WO Hz)^ 


1.2 X 10"^^ 



TABLE L Results of the cross-correlation of LIGO's HI and LI interferometers, analyzed for a potential power-law stochastic 
background of the form: f2gw(/) = ^^(//lOO Hz)™. The frequency range for each a is the band that contributes 99% of the 
full sensitivity, as determined by the inverse variance. All results correspond to the specified band, and an observation time 
of 218 hr. 90%-confidence Bayesian upper limits on Q,g^{f) (also expressed as limits on the strain noise density Sg{^{f)) are 
calculated from the point estimates and statistical errors, marginalising over a ±9% and ±15% uncertainty in the calibration 
magnitude of the HI and LI detectors. 



sensitive than Hl-Ll, due to the higher noise level of H2; 
the H2-L1 results are thus not used for the upper lim- 
its results (and are not shown), but they are consistent, 
within their error bars, with the Hl-Ll results. 

Systematic errors due to unresolved time variations in 
the interferometers' calibration and power spectra were 
investigated and determined to be small compared to the 
statistical error . Phase calibration uncertainties and 
timing errors are also negligible. Not negligible are inter- 
ferometer amplitude calibration uncertainties, estimated 
as a ±9% (±15%) frequency-independent uncertainty in 
the strain response magnitude for HI (LI). 

We construct a Bayesian posterior probability distri- 
bution for fla using the optimal point estimate fla and 
statistical error an^, j marginalising over the unknown cal- 
ibration magnitudes (see, e.g., [lOj). The prior probabil- 
ity distribution for is taken to be uniform from to 
0.02 (the maximum value corresponding to the largest 
background that is still consistent with the lowest single 
interferometer strain noise); the prior distributions for 
the calibration magnitudes are taken to be uniform be- 
tween 1 ± 0.09 (for HI) and 1 ± 0.15 (for LI). The 90% 
probability upper limit is then that value of Via for which 
90% of the posterior distribution lies between and the 
upper limit. (The upper limit is relatively insensitive to 
reasonable changes in the maximum value of the prior 
distribution for fio,, as well as to reasonable changes in 
the prior distributions for the calibration magnitudes.) 

The estimates for JIq, are entirely consistent with no 
stochastic background, within the sensitivity of the mea- 
surement. Furthermore, the cross-correlation spectrum 
(i.e., the integrand of Eq. ^ shows no distinct features 
(see Fig. and the « 26,000 values of Yj, with mean 
removed and normalized by the ayj , follow the expected 
normal distribution. 

Conclusions. -The energy density in a primordial 
background of GWs is constrained by big-bang nu- 
cleosynthesis theory, giving a conservative bound of: 
/ d(ln/)f]gw(/) < 1.1 X 10-5 0; if all the GW energy 
were concentrated in our sensitive band, this is still 60 x 
below the limit set here. A background from astrophys- 




50 100 150 

Cumulative analysis time (hr) 




100 110 120 130 
Frequency (Hz) 



140 150 160 



FIG. 3: The estimate Q^, for Hl-Ll and a = 0: (a) As a 
function of the amount of data analyzed. The shaded region 
shows the ±2(Tno band on Hq. (b) The real part of the Hl- 
Ll cross-correlation spectrum, in arbitrary units, (c) As a 
function of the frequency range analyzed. The shaded region 
shows the ±2(70p band on Q,q, cumulative in frequency from 
69 Hz to /. 



ical sources would be generated at much later cosmic 
times, and thus not subject to the above bound. In the 
LIGO band, such a background could be generated by 
the superposition of many short-lived signals, such as su- 
pernovae bursts or quasi-normal modes from compact ob- 
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jects (in the kHz region) and rotating neutron stars (in 
the tens to hundreds of Hz band). Uncertainties in the 
theoretical models are large, though the most optimistic 
predictions peak at i}g^{f) ~ 10"'^ 0- Nonetheless, the 
results presented here provide direct, measured limits to 
a stochastic background, that, in terms of energy den- 
sity, are nearly 10^ lower than previous measurements. 
Eventually, with 1-year of data at design sensitivity, the 
LIGO detectors will be sensitive at a level several times 
below the nucleosynthesis bound. 
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